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FIELDS 


D ECUS Program Library Write-up DECUS No. 8-191 

INTRODUCTION 


Ju Slds a demonstrat ion program which calculates and displays 
the surface potential of a given boundary conditioned plane. 

Each output facility is called by a 338 display push button giving 
a numeric and/or pictorial result. y y 


HARDWARE REQUIREMENTS 


PDP-8 with 8K store and 338 display 
HSR for Punch Display File 
Teletype for Print Matrix 


THEORY 


Given a point in a matrix it can be shown that the potential at 
that point approximates to the average of the potentials at the 
adjacent points. 


i.e. Vx,y= V*-l,y + Vx+l.y + Vx,y-1 + Vx,Y+l 
* 4 — 

The matrix is continually scanned with the above equation until the 
overall change is within the specified error. 

The calculation is more accurate if all eight adjacent points are 
used but the calculation time was found to be excessively long for 
a demonstration program. 

Reference for Theory: Modern Computing Methods NPL. Pages 117-118. 

Solution of Elliptical Differential Equations 

OPERATION 

Push Button Functions 

PB0 Redevelop display after each iteration 
PB1 Print results on Teletype 
PB2 Loop with display on 

PB9 Punch Display file on HSP 
PB10 Display angle (y increment 200) 

PBH Display angle (y increment 100) 

Operating Fields 


The Binary Tape contains the Floating Point Package eleminating the 
need for a separate load. 


Start Address = 200 









Starting the Push Button 2 set will cause the program to loop having 
initiated the display. Initially this is a simple vector formed:- 

Nil 

Resetting Push Button 2 causes the program to escape from the loop 
and enter the input routine. 


Data Input (via Teletype) 


This order of input must be strictly adhered to. 
Size of net - Example: 20 11 (integer) 


The program assumes an origin at 0,0 so the above input will give 
21 by 12 matrix points. The size of the net is limited by the 
storage available for the display file which starts in location 
01714. An input of 20 by 20 is close to the limit although the 
shape of the net will determine the real limit i.e. 30 by 10. 


Calculation error - Example: 0-01 (fractional) 

To read the error value the program enters the Floating Point Package. 


Boundary values - Example: 


20 0 301 
20 1 150.5 
20 2 10 
20 3 - 150.5 
20 4 -301 

The coordinate pairs are read by the programs own subroutine which 
will ignore LEADER, LF, and DELETE. 

An asterisk on the end of a numeric string will cause that string 
to be ignored, i.e. 29*20 0 301 = 20 0 301 

To read each boundary value the program enters the Floating Point 
Package. If a mistake is made in the value it can be corrected by 
overwriting, i.e. 


20 0 3011 
20 0 301 

The boundary point will be set to the last value read. 

Start Iteration - Example: 

This causes the program to escape from the input routine and begin 
the calculation. 

The push buttons should be set as required before giving this 
instruction. 
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OUTPUTS 


Display 


B e f°r e entering the calculation routine the program develops and 
initiates the display. This has a fixed width but the angle of 
view is dependent vertically on the setting of push buttons 10 and 11 . 

Example: Net size 5 3 


Push Buttons: 10, 11 



display routine to generate the final picture. See FIG.l. 

Pus ^ Button 0 is set, during the iteration, the display will be 
redeveloped after each pass. This enables the process of iteration 
to be viewed pictorially. 


Punch 


With Push Button 9 set the program enters the routine TAPE and 
punches out in binary format the complete display file. 

This can be read in with binary loader and displayed using the 
Push Button 3 Loop, allowing interesting results to be stored. 

Print 


This routine is entered when Push Button 1 is set and gives a print 
of the matrix on the Teletype. 

The output is broken into blocks of ten complete columns and con¬ 
tinues until each row is completed. 

The format can be altered by changing the output routine in the 
Floating Point Package. 

RESTART 


When the program has completed all requested operations it will re¬ 
turn to the input routine for further data. 

^ -*- s . squired to display the field at a different angle restart 
at 201 with push buttons 10 and 11 set as required. 
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GENERATING THE DISPLAY 

DISPLAY sets up the parallel line field as a series of vectors 
calculated from M,N and the push buttons 10 and 11. 

DCALC calculates the modified Delta-Y vectors in the format and 
does not effect the fixed X vectors. 
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INVISIBLE VECTORS 



)DYS)DYS) 
SHIFT VECTOR 
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DMAX “ DI1 


DELTA-Y Vector = Value (K) - Value (K-l) + DIJ 
When K = 0 Value (K-l) is forced to Zero 
When K = m Value (K) is forced to Zero 

This causes visible, vertical vectors to be generated at each end of 
the line when a + boundary condition is present. 
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CONCLUSION 


FIELDS demonstrates the graphical and numeric capabilities of the 
PDP-8 - 338 system. It has served as a useful introduction to the 
PDP-8 Floating Point Package and in simple 3-dimensional represent¬ 
ations of graphical results. 


Neither the results nor the display can be termed accurate and the 
scientist would much prefer an equipotential plot, a venture I have 
yet to consider. 
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APPENDIX Sample Data and Results 

20 20 


0.1 

00 20 0 

1 20 39 

2 20 77.2 

3 20 113 

4 20 147 

5 20 177 

6 20 202 

7 20 223 

8 20 238 

9 20 247 

10 20 250 

11 20 247 

12 20 238 

13 20 223 

14 20 202 

15 20 177 

16 20 147 

17 20 113 

18 20 77.2 

19 20 39 

20 20 0 
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20 0 0 


20 1 

-39 

20 2 

-77.2 

20 3 

-113 

20 4 

-147 

20 5 

-177 

^ 20 6 ■ 

-202 

20 7 ■ 

-223 

20 8 ■ 

-238 

20 9 - 

-247 

20 10 

-250 

20 11 

-247 

20 12 

-238 

20 13 

-223 

20 14 

-202 

# 20 15 

-177 

20 16 

-147 

20 17 

-113 

20 18 

-77.2 

20 19 

-39 


20 20 0 
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+ 
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192 
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CALCULATE 

DISPLAY 

VECTORS 


PRINT 

NET 

RESULTS 


(punch 

DISPLAY 
I FILE 


START 

AGAIN 
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